Detecting gravitational waves from inspiraling binaries with a network of 
detectors : coherent strategies by correlated detectors 
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We discuss the coherent search strategy to detect gravitational waves from inspiraling compact 
binaries by a network of correlated laser interferometric detectors. From the maximum likelihood 
ratio statistic, we obtain a coherent statistic which is slightly different from and generally better 
than what we obtained in our previous work. In the special case when the cross spectrum of two 
detectors normalized by the power spectrum density is constant, the new statistic agrees with the 
old one. The quantitative difference of the detection probability for a given false alarm rate is also 
evaluated in a simple case. 

PACS numbers: 95.85.Sz,04.80.Nn,07.05.Kf,95.55.Ym 



I. INTRODUCTION 

Gravitational waves (GW) is one of the most im- 
portant predictions of general relativity. In order 
to detect GW, several GW detectors are currently 
in operation around the world. However, the direct 
detection of GW has not been possible so far. Since 
the expected amplitude of GW is very small, multi- 
detector searches for GW are very important. A 
multi-detector search would (i) improve the detec- 
tion efficiency; (ii) improve our confidence in detec- 
tion of a GW event, if a candidate event is regis- 
tered; (iii) provide useful directional information on 
the GW event if the detectors are sufficiently geo- 
graphically separated; (iv) provide polarization in- 
formation if differently oriented. 

In our previous paper [lj (hereafter, Paper I), 
we discussed two different multi-detector detection 
strategies; the coherent strategy and the coincident 
strategy. The coherent strategy was discussed by 
Pai, Dhurandhar and Bose @, and by Finn Q. The 
target signal was the GW from inspiraling compact 
binaries for which the wave forms have been well 
studied. We considered the situation of two co- 
aligned and co-located detectors. We constructed 
the coherent statistic based on the maximum like- 
lihood method. We then compared the coherent 
method with the conventional coincident method 
by comparing the receiver operating characteristic 
(ROC) curves. We found that we could obtain better 
detection probability by the coherent method than 
by the coincident method. This conclusion was not 
changed even if the two detectors were correlated. 

In this paper, we consider two detectors with the 
same configuration as in Paper I. We discuss the re- 
fined derivation of the coherent statistic from the 



maximum likelihood statistic in the case when the 
detectors' noises are correlated. This statistic is 
slightly different from, and better than the one, we 
obtained in our previous work (Paper I). Here we 
present the mathematical proof of the superiority of 
the new statistic over the old one. In the special 
case when the cross-spectrum of two detectors nor- 
malized by the power spectrum density is constant, 
the new statistic agrees with the old one. Further, 
the quantitative difference of the detection probabil- 
ity for a given false alarm rate is also evaluated for 
a simple case. 



II. THE NEW COHERENT STATISTIC 

We consider two detectors which have same orien- 
tation, same arm length and the same noise power 
spectral density. We denote the one sided noise 
power spectrum density (PSD) common to each of 
the detectors as Sh(f) which has the property, 

(ni(f)nUf)) = ls h (f)S(f-f), (1) 

where nj(f) is the noise in detector / = 1,2. The 
angular brackets denote ensemble average. 

We assume the correlation of two detectors ex- 
pressed by 

(ni(f)n* 2 (f)) = ±e(f)S h (f)6(f-f), (2) 

where e(/)<SVi(/) represents the cross-spectrum. The 
function, e(/), is the cross-spectrum normalized by 
the power spectral density. We assume, for sim- 
plicity, the cross-correlation in the time domain is 
a function of \t — t'\. Thus, e(f) is real. 
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The response of a detector I to a gravitational 
wave from an inspiraling compact binary in the fre- 
quency domain is written as: 

Sj(f) =^ J (^^V;0/" 7/6 e i * (/;to ' 5o ' To ' T3) , (3) 

where t c and S c are respectively the coalescence time 
and the coalescence phase of the binary, tq and T3 are 
defined through mass parameters by Eq. (2.6) in Pa- 
per I. Ei((j), 9, l) is the extended antenna pattern 
function which depends on the orientation angles. 
The definitions of Ej and the amplitude normaliza- 
tion AT are given in Paper I. 

The function t c , S c , tq, 73) describes the phase 
evolution of the inspiral waveform. We adopt the 
3PN formula given by 

*(/; t e , S e , r , t 3 ) = 2irft c - 6 C - 1X - 
+ 12^ (7rM/) ~ 5/3 S a ^ M f) k/3 > ( 4 ) 

' fc=0 

where M is the total mass, rj is the ratio of the re- 
duced mass to the total mass, and ctk are given by 
Eq. (2.8) in Paper I. 

The form of ([3]) allows us to write the explicit 
quadrature representation of the i-th template as, 



Si(f;Pi,t c ,S c ) 



A(s (f; /2;,i c ) cos <5 C 

+ S7r/2(/;^i: f c)sin S c 



(5) 



where 



5o(/;Mi,«c) = a- 1 /- 7/6 e W(/;t -^ =0 ' To ' T3) , 



V2(/;/?i,*c) 
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/•/« r-7/3 



1/2 



(6) 
(7) 



Because the detectors are identically oriented, Ej = 
E. The templates so and have been normal- 
ized such that (so, so) = s^iz) — 1, where the 
scalar product (a, b) of two real functions a(t) and 
&(t) is defined as: 



•7* 



S*(/)6(/) 



&(/) 



(8) 



The quantities /j and /„ are respectively, the lower 
and upper cut off frequency. 

The output data from each detector, xi (I — 1, 2), 
can be written as xj = nj + sj, where ni is noise 
and sj is gravitational wave signal. We define two 
pseudo data by 



x± 

n± 
s± 



(xi ± x 2 ) 
n± + s±, 
(ni ± n 2 )/\/2, 
(ai±s 2 )/V2. 



(9) 
(10) 

(11) 



The two noise, n±, have property 
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(n ± (/)<(/')) = -(l±e(f))S h (f)S(f-f) 

= \s h(±) (J)5(f-f), (12) 



(rc + (/K_(/)) = «(/)n_(/))=0. 



(13) 



Since ni(f) are Gaussian variables, n±(f) are also 
Gaussian. Then, from Eq. (|13p . we find that n+(/) 
and ri-(f) are independent each other. Thus, it 
is straightforward to find the probability distribu- 
tion function of n±(t). We note that the signal, s_ 
is identically zero since the two detectors have the 
same orientation, the same arm length. We thus 
need to consider x + only. 

We derive the likelihood ratio from x + . We de- 
note the probability distribution function of x+(t) 
in the absence of signal as P(x+|0). The probability 
distribution function of 2+ (t) in the presence of the 
signal s+(t) is given by P(x + \s + ) = P(x + — s+|0). 
Since x+ is Gaussian in the absence of signal, the 
logarithm of the likelihood ratio, A + is given by, 



In A 1 



In 



P(x+\0) 



= (x+,s+) {+) - -(s+,s + ) (+) , (14) 

where the new inner product, (a, &)(+)> for real func- 
tions, a(t) and b(t) is given by, 

(a,b) (+) =2 df — . (15) 

Jf, b h (+)U ) 

The plus signal s + is given by, 
*+(/) = 

V2A (+) (s +) (/) cos S c + s { +l (/) sin <5 C ) , ( 16) 



where 



A+) =NEa {+)l 
rfu 



a 



(+) 



df 



f 



-7/3 



1/2 



(17) 
(18) 



and §n and 



f, Sh(+){f)_ 
^ l 2 are orthonormalized templates: 

(s +) ,4 +) ) (+) = (sW,sW) (+) =l, (19) 
(4 + \4/ 2 )(+)=0. (20) 
We define the normalization factor of s + as 

bf +) = (s + ,s + ) (+) =2Al +) . (21) 
The normalized s+ is given by, 

«+(/) = r^+Cf) 
fe (+) 

= 4 +) (/)cos<5 c + sW(/)sinJ c . (22) 
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Then the log-likelihood ratio becomes: 
lnA+ = b i+) {x + ,s + ) {+) 



5*» 



The maximization of In A+ with respect to 
S c (see (H) gives: 



(23) 
and 



max lnA^ 



1 



(+)\2 1 

M+)J> 



2 ^ncw- 



(24) 



We define 



„(+) 



(x + ,4 +) ) (+) , ^ = (^i^) (+)l (25) 



,(+) 



(+)- 



and we obtain 



Anew — (Cq ) + (c^y^) • 



(26) 



The mean and variance of c[ + ^ (i — 0, 7r/2) are given 

by 



\/2^4(+) cos<5 c , 



4 +) > = v^ (+)1 



„(+) 



» 2 



1. 



(27) 
(28) 
(29) 



Since c| + ^ obeys the Gaussian distribution, the PDF 
of A now in the absence of signal (Ar + \ = 0) is given 
by po(Ancw) = exp (-A ncw /2) /2. This gives the 
false alarm probability by one template as 



r^ltcmplatc 
FA 



dA now p (A 

new J 



cxp (-A new /2) . 



(30) 



This is the false alarm rate in the case when the pa- 
rameters of the signal are known. However, in gen- 
eral, we do not know, a priori, the parameters of the 
signal such as M, r], and t c . We thus need to compute 
A + for a bank of templates. This increases the false 
alarm rate. As in Paper I, we introduce the effec- 
tive number of statistically independent templates, 
iVjnd. The false alarm probability for the template 
bank thus takes the form, 

PFa = AUexp(-A* ew /2). (31) 

in the presence of signal is given 



The PDF of A new 

by 

Pi (A ncw ) = 



1 



exp 



X/n 



1 



(A„ 



2 

V2A {+) v/A 



(32) 



The detection probability becomes 

dAnewPl (A nC w) 



^ncw / A * \ 
DE l yi newJ 

dA 



1 



A- 



n ew 2 



1 



(An 



V2A/+) y/A 



(33) 



III. SUMMARY OF PAPER I 

In Paper I, we did not consider "±" data but used 
"±" correlations instead. The results of the earlier 
analysis are the following. 

The false alarm probability when the data is 
passed through A^ n d independent templates is given 
by Eq.(3.15) of Paper I: 



^(A old )=iVi„dexp(- 



A* 



id 



2(l + e )y ' 



(34) 



and the detection probability as given by Eq. (3.16) 
of Paper I is: 



Pde(Kh) 



where 



dA 



Id 



2(1 + co) 



exp 



(A old + 2A 2 ) 



2(1 + co) 



x J [A 



V2A 



old 



(l + eo)7 ' 



CO 



a -If: 



u ^(/)r 7/3 



(35) 



(36) 



IV. A MATHEMATICAL PROOF 

Now, we prove that for a given false alarm rate, 
the detection probability by A nC w is better than 
A id in general. First, we note that by setting 
t = A id/( 1 + e ), t* = A* ld /(1 + e ), and A = 
A/y/1 + cq, Eg. (f3~4"| and (|55|) are rewritten as 



PFA(f) = ^indCXpl - 



dt 



-exp 



-- (t + 2A\ 



xl a (Aq 



(37) 



We find that the formulae for the false alarm rate, 
(|3~lj) and |37|) are same. Thus, for a given false alarm 
rate, we have the same threshold. The functional 
forms of the detection probability, (|3"3"|) and (f3"T)) are 
the same. The difference comes from the amplitude 
of signal, -4.r+) and Aq in (|33|) and (|37| . 

If //(/) is a non-negative real function (a PSD for 
instance) we define 'average' of the function /i to be: 



(38) 



Then in terms of the average we may write, eo 
(</)>■ 

From the definition (|38|) . it follows that 



*(+) 



((! + </))■ 



(39) 
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From the Schwarz inequality, it immediately fol- 
lows that, 

<(i + £ (/))- 1 >((i + e (/))) > i, 

^±1(1 + e ) > 1. (40) 
er 

The Schwarz inequality Eq. (|4T)|) then implies 
^ -^o! the equality holding if and only if, e(f) is 
constant. This shows that the detection probability 
by A now is larger than that by A Q id. 




False alarm probability 

FIG. 1: The ROC curve (detection probability versus 
false alarm probability) of the old (dashed line) and new 
(solid line) statistic in the case when e(f) — 0.7 in the 
frequency 50Hz < / < 70Hz and zero otherwise. 



shown in FigfT] We see that the new statistic gives 
improved detection probability. For the false alarm 
rate of 1CT 10 , P$% = 0.814 and P£™=0.840. The 
difference is 0.026. 



VI. CONCLUSION 

We have discussed the coherent search strategy to 
detect gravitational waves from inspiraling compact 
binaries with a network of laser interferometric de- 
tectors. Based on the maximum likelihood method, 
we have derived a coherent statistic for detection 
of gravitational waves from inspiraling compact bi- 
naries by two aligned, correlated detectors. This 
statistic is slightly different from the one which we 
obtained in our previous work. We found that, in 
general, we can obtain a better detection probability 
for a given false alarm rate by the new statistic than 
by the old one. In the special case when the cross- 
spectrum normalized by the power spectrum density 
is constant, the new statistic agrees with the old one. 
The quantitative difference of the detection proba- 
bility for a given false alarm rate was evaluated in a 
simple case, and we found that there was a signifi- 
cant difference. In Paper I, we had concluded that 
the coherent detection strategy is better than the 
coincident strategy. Since we can have an increased 
detection probability by the new statistic, this con- 
clusion does not change even in the correlated case; 
on the contrary, the result is strengthened. 



AN EXAMPLE 
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We examine quantitatively the increase in the de- 
tection probability produced by the new statistic. 
We use LIGO I noise PSD given in Paper I. We take 
N\ n d — 6 x 10 6 which is approximately the value in 
Paper I. We set A — 7 '. We assume the following 
noise correlation model: 



e(f) 



0.7 (50Hz <f< 70Hz) 
otherwise 



(41) 



In this case, we have eo = 0.05665, A/*Jl + eo = 
6.809, and A (+) = 6.882. The ROC curves are 
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